function [matrix] = SampPIC(lx,rx,uy,ly,Msize)
matrix = zeros(Msize);
scale = Msize /64;
 [x, y] = Circle(Msize/2,Msize/2,25*scale,Msize);
 for i = 1:length(x)
    matrix(x(i),y(i)) = 128;
 end
 
 
 [x, y] = Circle(Msize/2,Msize/2,18*scale,Msize);
 
 for i = 1:length(x)
    matrix(x(i),y(i)) = 0;
 end
 

matrix(uy:ly,lx:rx) = 255;

function [x y] = Circle(centerX,centerY,R,Msize)
 gridx = meshgrid(1:Msize);
 gridy = meshgrid(1:Msize)';
 [x y] = find((gridx-centerX).^2 + (gridy-centerY).^2 < R.^2);





